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Abstract 


Orthogonal  frequency  division  multiplexing  ( OFDM)  is,  with  substantial  progress  in  digital 
signal  processing,  becoming  an  important  part  of  the  telecommunications  arena.  The  most 
appealing  feature  of  OFDM  is  the  simplicity  of  the  receiver  design  due  to  the  efficiency  with 
which  OFDM  can  cope  with  the  effects  of  frequency-selective  multipath  channels.  Here,  we 
address  the  problem  of  channel  estimation  for  OFDM  systems.  Exploiting  receive  antenna 
diversity,  a  second-order  statistics-based  (SOS)  blind  technique  is  proposed.  Our  method 
differs  from  the  existing  SOS-based  techniques  in  that  channel  estimation  is  carried  out 
using  the  frequency  domain  (i.e.,  post-FFT)  signals  whereas  existing  methods  use  the  time 
domain  signals  (i.e.,  pre-FFT).  In  the  proposed  method,  channel  identifiability  is  guar¬ 
anteed  regardless  of  the  channels  zeros  locations,  so  long  as  any  roots  common  to  all  the 
diversity  channels  are  on  the  unit  circle.  Only  short  data  records  are  required  to  achieve 
good  performance.  For  PSK  transmission,  this  method  enables  channel  estimation  even 
from  a  single  OFDM  symbol  at  high  SNR.  Further,  when  only  a  small  number  of  subcarri¬ 
ers  can  be  used  for  channel  estimation  because  of  computational  complexity,  we  determine 
the  optimal  set  of  subcarriers  in  terms  of  estimation  accuracy. 


keywords:  OFDM,  Multipath  channel,  Antenna  diversity.  Identification,  Estimation, 
Phase-Shift- Keying,  Quadrature- Amplitude-Modulation. 
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Chapter  1 

OFDM  Systems 


Orthogonal  frequency  division  multiplexing  (OFDM)  has  become  the  standard  of  choice 
for  wireless  LAN’s  such  as  HIPERLAN/2  and  IEEE  802.11a;  it  has  been  adopted  in 
Europe  for  Digital  Audio  Broadcasting  (DAB)  and  Digital  Video  Broadcasting  (DVB), 
MMAC  in  Japan,  and  fixed  wireless;  and  is  being  considered  for  several  IEEE  802.11 
and  802.16  standards,  including  wideband  Metropolitan  Area  Networks  (MAN)  [1].  The 
popularity  of  OFDM  stems  from  its  ability  to  transform  a  wideband  frequency  selective 
channel  to  a  set  of  parallel  flat-fading  narrowband  channels,  which  substantially  simpli¬ 
fies  the  channel  equalization  problem.  Because  of  the  time-frequency  granularity  that  it 
offers,  OFDM  appears  to  be  a  natural  solution  when  the  available  spectrum  is  not  con¬ 
tiguous,  for  overlay  systems,  and  to  cope  with  issues  such  as  narrowband  jamming.  In  the 
multi-user  context,  this  granularity  also  accommodates  variable  quality-of-service  (QoS) 
requirements  and  bursty  data. 

In  this  report,  we  review  the  basics  of  OFDM,  and  study  the  channel  estimation 
problem.  Clearly,  issues  such  as  timing  and  carrier  recovery  are  important,  but  will  not 
be  treated  in  this  chapter.  Channel  estimation  usually  consists  of  an  acquisition  phase 
followed  by  a  tracking  phase.  In  the  acquisition  phase,  initial  (often  coarse)  estimates 
of  the  channel  parameters  are  obtained.  Finer  estimates,  as  well  as  tracking  of  small 
time-variations,  of  these  parameters  are  acquired  during  the  tracking  phase. 

The  origins  of  OFDM  can  be  traced  back  to  a  1957  multi-carrier  HF  analog  modem  (the 
Kineplex),  but  it  was  only  with  digital  implementations  via  the  FFT  and  the  introduction 
of  the  cyclic  prefix,  both  by  Weinstein  and  Ebert  [2],  and  seminal  analyzes  by  Cimini  [3], 
that  OFDM  became  practical.  Several  recent  books  and  papers  on  OFDM  cover  different 
aspects  of  signal  and  system  design  [4,  5,  6,  7,  8]. 

Notation:  We  will  let  x(n,  m)  denote  the  m-th  element  of  the  vector  x{n).  Superscripts 

and  ^  will  denote  conjugate  transposition  and  transposition.  F  will  denote  the  FFT 
matrix  with  (A:,  m)th  entry  exp{—j2TTkm/M).  Tb  denotes  the  trace  operator. 
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1.1  Basics  of  OFDM 


Assume  that  the  information  bearing  symbols  are  to  be  transmitted  at  the  rate  of  R  sym¬ 
bols  per  second  over  a  multipath  propagation  channel.  The  duration  of  each  symbol  is 
therefore  Tg  =  1/R.  If  the  delay  spread^,  of  the  channel  is  larger  than  about  10%  of 
the  symbol  duration,  then  the  received  signal  may  suffer  from  significant  inter-symbol 
interference  (ISI),  which  can  drastically  increase  the  symbol-error-rate  (SER)^  unless 
counter-measures  are  undertaken.  Such  a  channel  is  said  to  be  dispersive  or  frequency- 
selective.  There  are  two  main  approaches  to  cope  with  such  channels.  The  first  approach 
is  to  use  a  single-carrier  system  with  an  equalizer  at  the  receiver  to  compensate  for  the 
ISI,  which  spans  [rmax/Tg]  symbols.  The  implementation  of  the  equalizer  may  become 
very  challenging  for  channels  with  large  delay  spreads,  and  at  higher  data  rates.  The  sec¬ 
ond  approach  is  based  on  multicarrier  modulation,  such  as  orthogonal  frequency  division 
multiplexing  (OFDM).  Here,  we  focus  on  the  latter  approach. 

The  operational  principle  of  an  OFDM  system  is  that  the  available  bandwidth  is 
divided  into  a  large  number  of  sub-bands,  over  each  of  which  the  wireless  channel  can  be 
considered  non-dispersive  or  fiat-fading.  The  original  data  stream  at  rate  R  is  split  into 
M  parallel  data  streams,  each  at  rate  R/M.  The  symbol  duration,  T,  for  these  parallel 
data  streams  is  therefore  increased  by  a  factor  of  M,  i.e.,  T  =  MTg.  Conceptually,  each  of 
the  data  streams  modulates  a  carrier  with  a  different  frequency  and  the  resulting  signals 
are  transmitted  simultaneously  (in  reality,  a  single  modulator  is  used,  as  we  discuss  later). 
Correspondingly,  the  receiver  consists  of  M  parallel  receiver  paths.  Due  to  the  increased 
symbol  duration,  the  ISI  over  each  channel  is  reduced  to  \Tmax/ (AfT^)]  symbols.  Thus,  an 
advantage  of  OFDM  is  that,  for  frequency-selective  fading  channels,  the  OFDM  symbols 
are  less  affected  by  channel  fades  than  are  single-carrier  transmitted  symbols.  This  is  due 
to  the  increased  symbol  duration  in  an  OFDM  system.  While  many  symbols  during  a 
channel  fade  might  be  lost  in  a  single-carrier  system,  the  sjunbols  of  an  OFDM  system 
can  still  be  correctly  detected  as  only  a  fraction  of  each  symbol  might  be  affected  by  the 
fade.  On  the  other  hand,  if  the  channel  is  time-selective,  i.e.,  the  channel  impulse  response 
varies  significantly  within  the  OFDM  symbol  period,  then  the  channel  matrix  is  no  longer 
Toeplitz,  and  conventional  OFDM  would  fail. 

Since  multicarrier  modulation  is  based  on  a  block  transmission  scheme,  measures  have 
to  be  taken  to  avoid  or  compensate  for  interblock  interference  (IBI),  which  contributes 
to  the  overall  ISI.  OFDM  systems  can  be  categorized  in  the  way  they  handle  IBI.  In  the 
most  popular  system,  a  guard  time  is  introduced  between  consecutive  OFDM  symbols 
as  a  cyclic  prefix,  (CP),  i.e.,  the  tail  end  of  the  OFDM  symbol  is  prefixed.  The  length 
of  the  cyclic  prefix  is  chosen  to  be  larger  than  the  expected  delay  spread;  after  proper 
time  synchronization,  the  receiver  discards  the  CP  and  thus  the  IBI  is  eliminated.  Time- 

^The  delay  spread  is  a  measure  of  the  delay  of  the  longest  path  (or  last  echo)  with  respect  to  that  of 
the  earliest  path  (or  first  arrival)  typically  the  root-mean-square  (RMS)  delay  spread  is  used. 

^The  symbol  error  rate  is  the  rate  of  errors  made  during  the  symbol  detection  process  at  the  receiver. 
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guarding  by  zero-padding  the  OFDM  symbols  has  also  been  proposed  in  [9,  10].  The 
issue  here  is  one  of  turning  the  transmitter  on  and  off  and  increased  receiver  complexity 
vs.  the  increased  SNR  and  decreased  SER.  Comparisons  between  cyclic-prefixing  and 
zero-padding  OFDM  systems  may  be  found  in  [11]. 

To  achieve  high  resilience  against  channel  dispersion,  a  large  number  of  subcarriers  is 
required.  However,  the  implementation  of  a  large  number  of  modulators  and  demodulators 
can  be  very  complex,  both  in  terms  of  the  physical  size  of  the  radio,  as  well  as  the  difficulty 
of  locking  in  multiple  oscillators.  This  complexity  can  be  significantly  reduced  by  digitally 
performing  the  modulation  and  demodulation  using  the  discrete  Fourier  transform  (DFT) 
and  its  inverse  (IDFT)  [2].  An  efficient  implementation  of  the  DFT  may  be  obtained  by 
any  of  the  available  fast  Fourier  transform  (FFT)  algorithms. 

The  choice  of  the  OFDM  parameters  is  a  tradeoff  between  various,  often  conflicting, 
requirements.  The  length  of  the  CP  is  dictated  by  the  delay  spread  of  the  channel. 
Introduction  of  the  CP  entails  a  reduction  in  rate  (or  wasted  bandwidth),  as  well  as 
a  SNR  loss;  to  minimize  these  inefficiencies,  the  number  of  subcarriers,  M,  should  be 
large.  However,  a  large  number  of  subcarriers  induces  a  high  implementation  complexity, 
increased  sensitivity  to  frequency  offset  and  phase  noise  (since  the  sub  carriers  get  closer 
to  each  other  as  M  increases)  and  an  increased  peak-to-average  power  ratio  (PAPR).  M 
is  dictated  by  concerns  regarding  practical  FFT  sizes  as  well  as  the  coherence  time  of 
the  channel.  We  will  not  address  the  issue  of  practical  choice  of  OFDM  parameters  here; 
we  refer  the  reader  to  [3,  6,  8].  In  this  chapter,  we  address  the  crucial  issue  of  CFO 
estimation. 

We  confine  our  attention  to  OFDM;  generalized  schemes  which  also  convert  frequency- 
selective  channels  into  a  bank  of  fiat-fading  channels  exist  [12];  see  also  [13].  In  [12],  the 
OFDM  scheme  is  shown  to  be  the  optimal  precoder  which  uses  a  cyclic  prefix.  In  a 
multiple-user  setting,  OFDMA  has  been  shown  to  be  the  optimal  scheme  in  the  sense  of 
maximizing  the  SNR  for  each  user  [14]. 

1.1.1  OFDM  Modulation 

OFDM  modulation  consists  of  M  (usually  a  power  of  2)  sub-carriers,  equi-spaced  at  a 
separation  of  A/  =  B/M,  where  B  is  the  total  system  bandwidth.  All  sub-carriers  are  mu¬ 
tually  orthogonal  over  a  time  interval  of  length  T  =  1/A/.  Each  sub-carrier  is  modulated 
independently  with  information-bearing  symbols  (this  does  not  preclude  coding  across  the 
sub-carriers).  Each  OFDM  block  is  preceded  by  a  CP  whose  duration  is  usually  longer 
than  the  delay  spread  of  the  propagation  channel,  so  that  IBI  can  be  eliminated  at  the 
receiver,  without  affecting  the  orthogonality  of  the  sub-carriers.  Practical  OFDM  systems 
are  not  fully  loaded  in  order  to  avoid  interference  between  adjacent  OFDM  systems:  some 
of  the  sub-carriers  at  the  edges  of  the  OFDM  block  are  not  modulated;  these  subcarri¬ 
ers  are  referred  to  as  virtual  subcarriers  (VSC).  The  number  of  these  VSC  is  dictated 
by  system  design  requirements  and  is,  in  general,  about  10%  of  M.  Some  subcarriers. 
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other  than  the  VSC,  may  also  be  deactivated.  For  example,  when  channel  state  infor¬ 
mation  (CSI)  is  available  to  the  transmitter,  subcarriers  experiencing  deep  fades  will  be 
left  unmodulated.  Further,  synchronization  preambles  are  often  made  by  nulling  a  large 
number  of  subcarriers.  Indeed,  a  preamble  consisting  of  a  repetition  of  two  identical  slots 
is  obtained  by  nulling  all  the  odd  subcarriers  [?].  Here,  deactivated  subcarriers  will  be 
referred  to  as  null-subcarriers  (NSC).  The  set  of  NSC  includes  the  VSC,  whose  placement 
and  number  are  imposed  by  system  design;  the  number  and  placement  of  the  remaining 
NSC  are  controlled  by  the  system  user,  and  could  vary  across  the  OFDM  symbols.  Let 
=  {0, ...,  M  —  1}  denote  the  entire  set  of  sub-carriers,  and  let  (resp.  denote  the 
subset  of  M  that  contains  the  Kn  (resp.  Zn)  modulated  (resp.  null)  subcarriers  during 
the  nth  OFDM  symbol  or  block. 


gi^Ttvn 


w(n) 


s(n) 


x{n) 


max  order  Lr 


Figure  1.1:  Discrete-time  complex  baseband  representation 

The  discrete-time  block  diagram  of  a  standard  OFDM  system  is  depicted  in  Fig. 
1.1.  The  vector  modulating  the  entire  set  of  subcarriers  during  the  nth  block  can  then 
be  expressed  as  s(n)  :=  V„SA;(n),  where  Sjc(n)  is  the  JQ-element  vector  of  symbols 
transmitted  on  the  activated  subcarriers,  and  V„  is  the  M  x  Kn  matrix  whose 
th  entry  is  one  if  the  ^th  symbol  is  transmitted  on  the  m-th  subcarrier  during  the  nth 
OFDM  block,  and  is  zero  otherwise.  Matrix  V„  is  a  full-rank  submatrix  of  an  (M  x  M) 
permutation  matrix.  We  assume  without  loss  of  generality  that  the  symbols  are  zero  mean 
and  have  unit  variance,  i.e.,  jF|s(n,  m)p  =  1.  The  (Mx  1)  data  block  s(n)  is  first  precoded 
by  the  IFFT  matrix  The  resulting  (M  x  1)  vector  u(n)  =  /?„F'^s(n)  is  called  the 
time-domain  block  vector,  or  the  time  domain  OFDM  symbol.  We  have  also  introduced 
a  normalization  parameter  /3„  :=  y/M/Kn  to  ensure  that  the  transmitted  power  is  kept 
constant  regardless  of  the  number  of  active  subcarriers.  Next,  a  CP  of  length  Lcp  is 
inserted  by  replicating  the  last  Lcp  elements  of  each  block  in  the  front.  The  redundant 
block  vector  can  be  expressed  as 


u{n)  =  [u{n,  M  —  Lcp),  u{n,  M  —  Lcp  +  1), ...,  u{n,  M  —  1),  u{n,  0), ...,  u{n,  M  —  1)]^  . 

The  P(=  M +Lcp)  samples  of  each  block  are  then  pulse  shaped,  upconverted  to  the  carrier 
frequency,  and  transmitted  sequentially  through  the  channel. 
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1.1.2  Demodulation 

We  model  the  frequency-selective  channel  as  an  FIR  filter  with  channel  impulse  response 
(CIR)  h  —  [ho,  where  L  is  the  channel  order.  In  practice,  the  system  is  usually 

designed  such  that  L  <  Lcp  <  M.  We  assume  that  the  CIR  is  time-invariant  over  AT  >  1 
consecutive  symbol  blocks,  but  could  vary  from  one  set  of  N  blocks  to  the  next. 

The  received  signal  is  downconverted  to  baseband  and  sampled  at  the  rate  of  P  sam¬ 
ples  per  extended  OFDM  symbol.  We  will  index  these  samples  by  [— Lcp,...,M  —  1]. 
We  will  assume  that  time  synchronization  has  been  achieved.  Discarding  the  samples 
i  —  —Lcp, ...,  —1  is  known  as  discarding  the  cyclic  prefix.  The  noise-free  received  signal 
corresponding  to  the  nth  OFDM  symbol,  u(n),  can  be  written  as 

L 

x{n,  k)  —  ^2  heu{n,  k  —  £)  , 
e=o 

for  k  =  0, ...,  M  —  1.  Recall  that  with  the  insertion  of  CP,  we  have  u{n,  £)  =  u{n,  M  +  £) 
for  t  =  —Lcp, ...,  —1,  and  u{n,  i)  =  u{n,  £)  ior  £  =  0, M  —  1.  Then,  collecting  samples, 
A:  =  0, ...,  M  —  1,  of  x(n,  k),  we  obtain 


'ho  0  ...  0  0  ...  0  ’ 

x(n,0) 

hi  ho  ...  0  0  ...  0 

u{n,  0) 
u{n,  1) 

x{n,  M  —  1) 

hL-i  hi-2  —  •  ^0  0  •  •  •  0 

hi,  hi^i  hi  ho  0  0 

_  u{n,  M  —  1)  _ 

hi  /i2 

.  .  . 

hL-i 

hL  0  ••• 

0  ' 

1 - 

1  1 

^  ^  •• 

+ 

/l2  hs 

hi 

0  0  ••• 

0 

o  o 

^  O 

0 

0 

0 

0 

o  o 

o  o 

0 

0  . 

^  u{n,  M  -  Lcp)  \ 

ho 

0 

...  0 

0 

0  hi 

hi-i  . . . 

hi 

“ 

hi 

ho 

...  0 

0 

0  0 

hi 

/l2 

u(n,  0) 

* 

: 

i 

* .  : 

: 

I 

u{n,  1) 

hi-i 

Ai-2 

...  ho 

0 

0  0 

0  ... 

hi 

: 

hi 

hc-1 

...  hi 

ho 

0  0 

0  ... 

0 

_  u{n,  M  —  1)  _ 

x{n)  = 

Hcw(n) 

J 

Notice  that  the  matrix  He  is  circulant  wdth  first  column,  [ho,  hi,  ...hL,0,  ...,0]^. 

Under  the  narrowband  assumption,  the  noise  corrupted  counterpart  of  (1.1)  becomes 
(assuming  no  carrier  offset  i.e.  z/  =  0  in  Fig.  1) 

x(n)  =  /?„HcF^s(n)  -f-  n(n)  (1.2) 
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where  He  is  the  (M  x  M)  circulant  matrix  defined  in  (1.1),  and  n{n)  is  the  M  x  1  noise 
vector,  which  is  assumed  to  be  zero  mean  circular  Gaussian  with  covariance  matrix  cr^I. 

Demodulation  is  based  on  the  well-known  property  that  any  circulant  matrix  can  be 
diagonalized  by  pre-multiplication  by  the  FFT  matrix  and  post-multiplication  by  the 
IFFT  matrix.  Let  Hk  ,  denote  the  frequency  response  of  the  channel 

at  frequency  2'Kk/M,  and  let  H  :=  diag  (LTo?  •••,  fifM-i).  Then,  the  signal  in  (2.1)  can  be 
rewritten  as 


x{n)  =  /?„F"HFF^s(n)  +  n(n) 

=  pnF^Hs{n)  +  n{n)  .  (1.3) 

Therefore,  after  FFT  processing,  the  so-called  frequency-domain  received  symbol  blocks 
are  obtained  as 


y{n)  =  Fa3(n)  = /?„FF^Hs(n)  +  T7(n)  (1.4) 

where  rj{n)  =  Fn(n)  is  again  AWGN  with  covariance  matrix  a‘^1.  In  the  absence  of  carrier 
offset,  the  frequency-domain  blocks  are  then  obtained  as 

y{n)  =  ;0nHs(n)  -I-  'n{n)  .  (1.5) 

The  effect  of  the  frequency-selective  channel  on  the  OFDM  signal  is  completely  captured 
by  scalar  multiplications  of  the  data  symbols  by  the  frequency  responses  of  the  channel  at 
the  subcarrier  frequencies.  Further,  demodulation  at  the  receiver  does  not  color  the  addi¬ 
tive  noise.  If  none  of  the  channel  zeros  coincides  with  an  activated  subcarrier,  maximum 
likelihood  detection  of  the  symbols  is  straightforward.  Zero-forcing  and  MMSE  equalizers 
can  be  applied  on  a  per-carrier  basis.  From  eq.  (1.5),  we  see  that  under  the  constraint  of 
constant  transmitted  power,  the  presence  of  NSC  (i.e.,  =  M/Kn  >  1)  implies  a  higher 

local  SNR®  at  the  modulated  subcarriers  at  the  expense  of  bandwidth  efficiency.  But,  for 
frequency-selective  channels,  the  local  SNR  can  vary  significantly  across  the  subcarriers. 
Information  transmitted  on  a  subcarrier  that  is  experiencing  a  deep  fade  (i.e.,  low  SNR) 
could  be  lost,  i.e.,  frequency  selectivity  of  the  channel  could  degrade  BER.  There  are  two 
major  techniques  to  mitigate  this  problem.  The  first  approach  is  to  code  across  the  sub¬ 
carriers,  typically  by  using  a  convolutional  code.  This  improves  SER  at  the  expense  of 
reduced  rate,  but  does  not  require  CSI  at  the  transmitter.  The  second  approach  is  based 
on  power  or  bit  loading  techniques,  and  assumes  that  CSI  is  available  at  the  transmitter. 
System  capacity  can  be  maximized  by  adapting  the  powers  or  the  bit  loads  of  the  different 
subcarriers  to  the  channel.  In  the  case  of  bit  loading,  the  constellation  sizes  of  the  symbols 
transmitted  on  the  different  subcarriers  could  be  adjusted  according  to  the  corresponding 
SNRs  using  a  water-filling  method  [15].  This  is  the  t3q)ical  scenario  in  discrete  multitone 
(DMT),  as  the  wired  version  of  OFDM  is  called  [16]. 

®The  local  SNR  at  subcarrier  k  is  defined  as  {M/Kn){\Hk\^ /a^) 
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1.2  Channel  Estimation 


As  mentioned  previously,  the  main  advantage  of  OFDM  systems  is  the  reduced  complexity 
of  the  equalizer  at  the  receiver,  which  allows  for  inexpensive  hardware  implementation. 
Since  OFDM  transforms  a  frequency-selective  cha.nnel  into  parallel  flat-fading  subchan¬ 
nels,  a  bank  of  one-tap  equalizers  suffices.  Clearly,  the  channel  needs  to  be  estimated  in 
order  to  design  the  equalizer.  Channel  identiflcation  may  also  be  used  for  shortening  the 
channel  impulse  response  and  determining  power  loading  at  the  transmitter  (if  channel 
state  information  is  available  at  the  transmitter,  channel  nulls  can  be  avoided  by  not 
transmitting  symbols  on  those  subcarriers). 

Traditionally,  channel  estimation  is  carried  out  using  pilot  symbols.  Because  they  save 
bandwidth  and  are  capable  of  tracking  (slow)  channel  variations,  blind  channel  estima¬ 
tion  and  equalization  methods  are  well  motivated.  Numerous  blind  channel  identification 
algorithms  for  OFDM  with  redundancy^  have  recently  been  developed.  Some  of  these 
methods  are  based  on  the  cyclostationarity  of  oversampled  OFDM  signals  [28],  while  oth¬ 
ers  are  based  on  subspace  decompositions  that  exploit  either  the  redundancy  introduced 
by  the  cyclic  prefix  (CP)  [29],  or  require  changes  to  the  current  systems  such  as  the  use 
of  precoding  [30],  or  zero-padding  [31].  Some  of  these  methods  require  many  OFDM 
symbols  for  channel  estimation,  so  that  they  have  to  be  modified  for  use  over  rapidly 
fading  channels.  Moreover,  in  subspace  methods  for  CP-OFDM,  channel  identifiability 
is  not  guaranteed  if  the  channel  has  nulls  on  the  subcarriers  [29].  A  frequency  domain 
method  based  on  the  finite-alphabet  of  the  information-bearing  symbols  was  developed 
in  [32].  Blind  channel  estimation  of  single-input-multiple-output  (SIMO)-OFDM  models 
has  apparently  not  been  studied.  In  the  next  chapter,  we  exploit  receive  antenna  diversity 
and  develop  a  blind  identification  technique  using  the  frequency-domain  signals. 


^Blind  channel  estimation  for  OFDM  without  redundancy  has  been  addressed  in  [27]. 
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Chapter  2 


Blind  Channel  Estimation  using 
Receive  Space  Diversity 

2.1  Introduction 

The  use  of  multiple  receive  antennas  has  been  largely  in  the  context  of  diversity  combining, 
e.g.,  [33,  34],  where  the  channel  is  assumed  known,  or  is  estimated  via  pilot  symbols,  or 
spatial  diversity  is  used  only  for  beamforming  and/or  channel  tracking.  Blind  channel 
estimation  of  SIMO-OFDM  models  has  apparently  not  been  studied.  In  this  report,  we 
exploit  receive  antenna  diversity  and  develop  a  blind  identification  technique  using  the 
frequency-domain  signals.  Our  method  possesses  the  following  attractive  properties: 

•  It  is  performed  in  the  frequency  domain.  Therefore,  unlike  the  time-domain  meth¬ 
ods,  one  can  control  the  number  and  frequencies  of  the  subcarriers  to  be  used  in 
channel  estimation.  We  will  later  show  how  to  select,  for  a  fixed  number  of  subcar¬ 
riers,  the  optimal  set  of  subcarriers  in  terms  of  estimation  accuracy.  On  the  other 
hand,  properly  designed  pre-FFT  equalizers  could  reduce  the  number  of  required 
FFT  chipsets. 

•  Unlike  the  existing  methods,  it  does  not  suffer  from  the  self  noise  effect  when  the 
symbol  constellations  are  constant  modulus. 

•  It  is  SOS-based  unlike  the  method  in  [32]  which  uses  higher-order  statistics.  For  ex¬ 
ample,  for  MPSK  symbols,  the  method  in  [32]  requires  M*'’-order  statistics.  Its  per¬ 
formance  deteriorates  with  M  whereas  that  of  the  proposed  method  is  independent 
of  M.  Further,  our  method  can  cope  with  a  larger  delay  spread,  i.e.,  L  <  {K  —  l)/2 
where  K  is  the  number  of  modulated  subcarriers,  compared  with  the  method  in  [32] 
for  which  L  <  (K  —  1)/M.  The  latter  condition  might  not  be  satisfied  for  large 
constellation  sizes. 

•  It  is  robust  to  channel  order  over-determination. 
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2.2  Signal  Model 

In  this  chapter,  the  ^th  coefficients  for  the  zth  channel  is  denoted  by  hi{C).  After  dis¬ 
carding  the  cyclic  prefix,  FFT  processing,  and  ignoring  the  null-subcarriers  (NSC’s),  the 
frequency-domain  nth  OFDM  block  at  the  zth  receive  antenna  can  be,  using  the  signal 
model  presented  in  Chapter  1,  written  as 

V^{n)  =  B.iS{n)-\■'ni{n),  n  =  i  =  l,...,/,  (2.1) 

where  y(n)  is  the  [K  x  1)  vector  whose  entries  are  y{n,  k),  k  e  )C,  r}{n)  is  additive  noise, 
with  Hi  =  diag  (i/i,fc,  k  e  JC),  and 

H,,,  =  '£hi{e)p;‘ ,  p*  :=  . 

e=o 

In  the  above  model,  we  have  assumed  the  /9„’s  defined  in  Chapter  one  to  be  identical 
across  the  blocks  and  they  are  normalized  (with  respect  to  the  noise  power)  to  unity 
Notation.  We  will  let  superscript  *,  and  1  denote  conjugate,  conjugate  transpose, 
and  pseudo-inverse;  8mn  =  1  if  m  =  n  and  0  otherwise;  5ijki  =  1  if  i=j=k=i  and  zero 
otherwise;  Ik  is  the  K  x  K  identity  matrix. 

We  will  assume  that  i)  the  input  sequence  s(n)  is  zero  mean  and  E  {s(n)s-^(m)}  = 
SnmJ-K,  and  m)  the  additive  noise  is  uncorrelated  with  the  data  sequence,  spatially  and 
temporally  white  and  circularly  symmetric  Gaussian  (STCWG)  with  variance 

2.3  Identifiability  Conditions 

The  cross-correlation  between  yi{n,k)  and  yj{n,k),  i.e.,  between  the  A:-th  subcarriers  at 
antennas  i  and  j  is  given  by 

'tijik)  :=  E  {y,(n,  k)y*{n,  A;)}  =  -f  (2.2) 

with  k  E  JC,  and  i,j  =  1, ...,  /,  and  where  we  have  used  the  fact  that  E {|s(n,  A:)p  =  1}. 


2.3.1  Identifiability  of  the  Additive  Noise  Power 

If  the  additive  noise  power  is  unknown,  it  can  be  estimated  by  exploiting  receive  diversity. 
Two  diversity  branches  are  sufficient  to  identify  cr^.  Indeed,  using  jij{k),  i  ^  j,  we  readily 
obtain  the  following  equation 

-  hi{k)  +  Jjj{k)]a^  -t-  7ii{k)jjj{k)  -  |7ij(^)P  =  0  . 

Thus,  cr^  can  be  obtained  as  the  unique  solution  to  a  quadratic  (only  one  root  is  valid). 
Combining  all  the  correlations,  we  can  obtain 


-T 


a{k)  —  ^a{k)^  —  4/3(A:) 


fceAC 
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where  a{k)  =  _  ^  buik) +'rjj{k)] 


and  P(k)  = 


A.. 


Y1  bii{k)ljj{k)  -  bijik)]"^ 


Note  that  two  diversity  branches  are  sufficient,  and  that  cr^  can  be  identified  regardless 
of  the  values  of  the  channel  coefficients,  the  i?i(pfe)’s.  Here,  for  simplicity,  we  focus  on 
the  case  of  temporally  and  spatially  uncorrelated  noise  keeping  in  mind  that  most  of  the 
results  established  here  can  be  extended  to  the  case  of  temporally  correlated  but  spatially 
uncorrelated  noise. 


2.3.2  Channel  Identifiability 

The  first  term  of  the  RHS  of  eq.  (2.2)  can  be  expressed,  using  the  z-transform,  as 

where  Gij{z)  =  Hi{z)Hj{l/z*)  and  Hi{z)  =  Y^e^ohi{£)z~^.  The  inverse  z-transform  of 
Rij{z)  =  Gij{z)  -I-  cr^  is  given  by 

min{L,L—r) 

'^ijb)=  X)  hi{l  +  T)h*(l)  +  a'^5ij5{T),  (2.3) 

;=m£ix(0,-r) 

with  r  =  —L, ...,  L. 

Let  Hij  be  the  Toeplitz  matrix  with  {m,n)  entry  rij{m  —  n);  then, 

Rij  =  niHf  +  (2.4) 

where  Hi  is  the  (L-f-lx2I/-fl)  Toeplitz  matrix  with  first  row 

[/ii(0), ..., /ii(L),  0, ...,  0]^  and  first  column  [/ij(0),0,  ...,0]^.  Concatenating  the  Rjj’s,  we 
obtain  the  following  I{L  +  1)  x  I{L  +  1)  matrix  R  =  which  can  be  expressed  as 

R  =  HH^  +  (2.5) 

where  H  =  [Hi...Hjb.  Notice  that  R  is  the  same  as  the  theoretical  covariance  matrix  of 
a  single  carrier  SIMO  system.  This  implies  that  all  the  existing  SOS-based  channel  iden¬ 
tification  methods  could  be  applied  to  our  matrix  R  provided  the  latter  can  be  uniquely 
estimated  from  the  jij{kys.  This  condition  is  satisfied  if  the  number  of  modulated  sub¬ 
carriers  K  >  2L. 

If  R  is  uniquely  identified,  then  the  issue  of  channel  identifiability  is  similar  to  that  of 
single  carrier  SIMO  systems.  The  conditions  for  data  recovery^  are  however  different  from 

^By  data  recovery  we  mean  that  in  the  absence  of  noise,  the  data  streams  can  be  exactly  determined 
up  to  a  sign  ambiguity. 
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those  of  single  carrier  SIMO  systems  (without  cyclic  prefix).^  In  the  literature,  channel 
identifiability  is  often  amalgamated  with  data  recovery.  Here,  we  distinguish  between  the 
two  concepts.  The  following  Proposition  states  the  channel  identifiability  conditions  when 
SOS  are  used. 

Proposition  1.  The  {L  +  l)-tap  channels  of  an  I  receive-diversity  OFDM  system,  in 
the  presence  of  STCWG  noise,  can  he  blindly  identified  using  second-order  statistics  (up 
to  a  phase  ambiguity)  if  the  number  of  modulated  subcarriers  K  >  2L,  and  any  roots  that 
are  common  to  all  the  channels  are  on  the  unit  circle. 

Proof.  Once  the  i2jj(z)’s  and  are  identified,  the  Gij(z)’s  can  be  determined  as 
Gij(z)  =  Rij{z)  —  o^.  If  0  is  a  root  of  Gii{z)  then  its  reciprocal  1/6*  is  also  a  root 
of  Gii{z).  This  implies  that  identifying  Gii(z)  leads  to  spectrally  equivalent 

channels  Hi{z),  where  q{i)  is  the  total  number  of  roots  Hi{z)  ^  and  qi{i)  is  the  number  of 
unit-modulus  roots  of  Hi{z).  The  ambiguity  in  the  non  unit-modulus  (NUM)  roots  can 
be  eliminated  using  receive  diversity  as  shown  next.  Note  that  at  least  one  of  the  /ii(0)’s 
has  to  be  non-zero  to  account  for  the  inherent  delay  ambiguity  in  blind  techniques. 

For  simplicity,  we  assume  that  hi{0),hj(0)  7^  0^.  If  and  denote  the  roots  of 
Hi{z)  and  Hj{z)  respectively,  then 

?(*)  9(j) 

G«(z) = h,(o)h-(o)  nil  -  Hi'  -  0^1 

1=1  1=1 

Therefore  identifying  the  NUM  roots  of  Gij(z)  will  provide  the  set  of  roots 
I  =  l,...,5i(i);  m  =  l,...,qi{j).  If  Hi{z)  and  Hj{z)  share  no  common  NUM  roots,  then 
the  channels  Hi{z)  and  Hj{z)  can  be  blindly  identified  (up  to  a  phase  ambiguity)  using 
Gij{z),  Gij{z)  and  Gjj{z).  More  generally,  if  no  NUM  root  is  common  to  all  the  channels, 
then  all  the  channels  can  be  uniquely  identified  (up  to  a  phase  ambiguity)  using  iJy, 
i,j  =  1,...,/.  These  results  are  also  valid  for  single  carrier  SIMO  systems  provided  the 
covariance  matrix  is  uniquely  identified.  □ 

Note  that  for  single-carrier  systems,  channel  identifiability  is  guaranteed  under  the 
same  conditions  as  above  except  for  the  condition  K  >  2L  which  is  not  applicable. 

2.3.3  Data  Recovery 

For  single  carrier  systems,  data  recovery  is  guaranteed  if  the  matrix  TL  is  full  column  rank, 
which  is  obtained  if  the  channels  share  no  common  zeros.  For  our  OFDM  system,  data 
recovery  is  guaranteed  under  a  weaker  condition  as  stated  next. 

Proposition  2.  For  an  I  receive-diversity  OFDM  system  operating  over  a  {L-\-l)-tap 
channel,  in  STCWG  noise,  data  recovery  is  guaranteed  using  second-order  statistics  (up 

^Data  recovery  conditions  and  channel  estimation  for  cyclic-prefixed  single  carrier  systems  are  similar 
to  those  of  our  OFDM  system. 

®if  hi{l),hi{L)  ^  0  then  g(i)  =  L 

“^We  can  check  whether  hj{0)  7^  0  by  inspecting  the  ry’s  in  eq.(2.3). 
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to  a  phase  ambiguity)  if  the  number  of  modulated  subcarriers  K  >  2L,  any  roots  that  are 
common  to  all  the  channels  are  on  the  unit  circle,  and  the  common  unit-modulus  roots,  if 
any,  do  not  coincide  with  the  modulated  subcarrier  frequencies,  2'Kk/N,  k  E  JC. 

Proof.  The  first  conditions  are  required  for  channel  identifiability.  Once  the  channels 
have  been  identified,  their  frequency  responses  at  the  subcarriers  frequencies  are  calcu¬ 
lated.  If  the  channels  have  a  common  root  that  is  constant  modulus  and  its  phase  coincides 
with  one  of  the  subcarriers  frequencies,  then  the  corresponding  frequency  responses  will 
be  zero  at  this  subcarrier,  in  which  case  data  cannot  be  recovered.®  □ 

2.4  Channel  Estimation  Approach 

Channel  estimation  is  based  on  the  estimation  of  the  matrix  R,  which  can  be  done  either 
in  the  time-domain,  or  based  on  the  frequency  domain  correlations,  the  7y’s.  Here,  we 
use  the  latter  approach. 

We  first  estimate  the  7^j(A:)’s  for  A:  €  /C  and  i,j  =  1, ...,/, 

1  ^ 

%  ;  *>•?  =  1’  (2.6) 

n=l 

where  ©  denotes  the  Shur-Hadamard  product  (i.e.,  element  wise  multiplication).  Esti¬ 
mates  of  the  Vij  =  [rij{—L),  ...,rij{L)]'^  are  obtained  via 

=  =  -.I  (2.T) 

where  Bx;  is  the  K  x  {2L  1)  matrix  whose  elements  are 

exp{j2Trkl/N)  with  k  E  JC,  the  set  of  active  subcarriers,  and  I  =  —L, ...,  L. 

For  constant-modulus  symbols  and  in  the  absence  of  noise,  %j{k)  =  'yij{k)  so  that 
Ry  =  Rjj.  This  implies  that  for  PSK  constellations,  the  channels  can  be  exactly  identified 
(up  to  a  phase  ambiguity)  using  a  single  OFDM  block  in  the  absence  of  noise.  Channel 
identification  using  the  time  domain-based  estimation  of  covariance  matrix  R  does  not 
benefit  from  this  nice  property.  It  is  this  property  which  makes  the  proposed  method 
particulary  efficient  for  PSK-OFDM.  (In  the  absence  of  noise,  it  was  shown  in  [35]  that 
both  channel  and  S3mabols  can  be  recovered  from  a  single  OFDM  symbol  with  m-PSK 
modulation  if  #  sub-carriers  N  >  mL.)  For  non-constant  modulus  symbols,  e.g.,  QAM 
constellations,  M  is  required  to  be  moderate/large  in  order  for  7y  (A;)  to  be  a  good  estimate 
of  jij(k)  even  in  the  absence  of  noise.  The  is  referred  to  as  the  self  noise  effect. 


2.5  Performance  Analysis 

Here,  we  derive  the  covariances  of  the  cross-correlation  estimates.  Once  these  covariances 
are  obtained,  performance  analysis  of  channel  estimates  can  be  carried  out  along  the  lines 

®Note  however  that  coding  across  the  subcarriers  is  always  used  in  order  to  address  this  problem. 
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of  [36].  Because  the  estimation  of  the  cross-correlations  is  carried  out  in  the  frequency- 
domain,  we  are  able  to  provide  closed-form  expressions  for  the  finite-sample  covariances 
of  the  cross-correlations. 

5.1  Channel-dependent  bounds 

The  finite-sample  covariances  of  the  sample  cross-correlation  7^ based  on  M  OFDM 
blocks,  are  given  by 

Mcov  {jij ,  Jk,e)  =  Hf  -I-  a^Sje^iRk 

a^6ikRfJle  +  (^ijkeCin  +  Sik^je)o’^^K 

where  pis  =  E  {\s{n,rn)\^}  /{E  {\s{n,Tn)\^})^,  and  c^r,  ■=  —2  -\-E{\r}\‘^}/a^  is  the  (nor¬ 
malized)  fourth-order  cumulant.  Note  that  cov  {'yij,'y*k^e)  =  cov  (7ij,7£_fc)  since  7fc,^(.)  = 

For  constant-amplitude  symbols,  e.g.,  PSK  constellations, 

P4s  =  l,  and  the  first  term  on  the  RHS  of  the  above  equation  vanishes.  This  implies  that 
for  PSK  symbols,  the  cross-correlations  estimates  are  asymptotically  (i.e.,  as  cr^  — >•  0) 
consistent.  This  is  in  contrast  with  the  case  of  QAM  constellations  for  which  p4s  ^  1. 
The  corresponding  cross-correlation  estimates  are  consistent  only  when  the  number  of 
blocks  M  tends  to  00.  Note  also  that  {p4s  —  1)  increases  with  the  constellation  size,  and 
so  do  the  errors  on  the  cross-correlation  estimates.  For  instance,  p4s  =  1.32  for  16QAM 
and  p4s  =  1.38  for  64QAM  (with  a  limiting  value  of  1.4).  Using  eq.  (2.7),  the  covariances 
of  the  sample  estimates  of  are  obtained  as 

cov  {r-ij,  rlk)  =  cov  (f^ ,  rk,i)  =  Bj[;Cov  (tij,  7k.«)  (BJ;)®'  •  (2.8) 

5.2  Channel-independent  bounds 

These  covariance  expressions  are  channel  dependent.  As  wireless  channels  are  random, 
it  is  interesting  to  average  the  covariance  expressions  over  all  the  possible  channels.  We 
assume  that  the  channels  are  independent  and  Rayleigh  distributed  (the  Ricean  case 
could  be  treated  in  the  same  manner).  Let  Q  denote  the  covariance  matrix  of  the  channel 
coefficients,  the  hj’s.  For  simplicity,  we  focus  on  the  case  of  uncorrelated  scattering,  for 
which  Q  is  diagonal;  then,  with  ajj  =  Tr  {Q},  we  obtain 

E  . 

In  the  case  of  uncorrelated  scattering,  the  averaged  finite-sample  covariances  of  the 
sample  cross-correlation  j  are 

co\/ {‘yij,^k,e)  =ocijkelK,  where 

^ijkt  ~  ^ ^\^ij^kt  4"  ^ik^i(\  4"  2(Tjj<T  S^k^ji 

4“  <7  SikSjg  "h  Sijk£C4'q  . 
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The  averaged  errors  on  the  estimation  of  the  are  given  by 

COV(fy,ffc^)  =  . 


5.3  Optimal  Selection  of  Subcarriers 

In  order  to  reduce  the  computational  complexity  of  the  proposed  channel  identification, 
it  is  desirable  to  reduce  the  number  of  subcarriers  used  for  the  estimation  of  the  cross¬ 
correlation  coefficients.  In  this  subsection,  we  answer  the  following  question:  for  a  given 
number  of  subcarriers  to  be  used  in  bhndly  estimating  the  channels,  what  is  the  optimal 
set  of  subcarriers  in  terms  of  estimation  accuracy. 

In  the  previous  subsection,  we  have  shown  that  under  the  uncorrelated  scattering 
scenario,  the  errors  in  the  estimation  of  the  7^-  are  independent  of  the  frequencies  of  the 
subcarriers  set.  The  errors  in  the  estimation  of  the  do  however  depend  on  K. 

Proposition  3:  The  optimal  selection  of  D  subcarriers  in  terms  of  channel  estimation 
performance  is  given  by 

V  =  arg  min  TV 

VCK. 

Thus,  if  D  <  K/2,  then  the  optimal  set  V  is  obtained  by  selecting  equally  spaced  subcar¬ 
riers  with  maximum  spacing  between  the  subcarriers.  □ 

2.6  Channel  Estimation  Algorithms 

Given  estimates  of  fij,  the  noise  power  can  be  estimated  using  results  in  Sec  3.1. 

Cross-Correlation  Matching  Method:  Let  fij  =  —  cr^^y-ei+i  and  fij  —  Vij  — 

where  e£,+i  is  a  (2L  -f-  l)-element  vector  whose  (L  -I-  l)th  element  is  one,  and 
the  rest  are  zeros.  Now,  we  propose  to  estimate  the  channels  using  the  following  matching 
criterion 

I 

J(h)=  Y.  (2“) 

where  Wy  are  weighting  matrices  and  |e|wij  =  e^Wije  for  any  vector  e  of  appropri¬ 
ate  dimensions.  The  optimal  weighting  matrices  can  be  obtained  using  eq.  (2.8).  The 
performance  of  COM  can  be  derived  along  the  lines  of  [36],  which  will  then  serve  as  a 
benchmark  for  other  SOS-based  methods. 

Subspace  Method:  Using  eq.  (2.5),  we  can  readily  adapt  the  subspace  approach  of 
[37]  to  obtain  channel  estimates.  Note  the  slightly  relaxed  conditions  of  Propositions  1 
and  2  in  the  context  of  multi-carrier  signals.  Note  that  the  deterministic  subspace  method 
does  not  suffer  from  self  noise  even  for  QAM  constellations.  However,  the  identifiability 
condition  for  these  methods  is  that  the  channels  have  no  common  root,  which  is  more 
restrictive  than  the  identifiability  condition  established  above.  Moreover,  Deterministic 
subspace  methods  fail  if  the  channel  order  is  not  estimated  properly. 
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Root  Selection  Method:  The  proof  of  Proposition  1  outlines  the  method.  Prom 
Gij  we  obtain  the  roots  of  Hi{z)H2{l/z*).  Every  root  of  G12  can  be  compared  with  those 
of  Gii,  and  the  closest  selected.  The  set  of  L  roots  of  Gn  closest  to  those  of  Gu  yields 
Hi{z).  H2{z)  is  obtained  similarly. 

2.7  Simulations 

Since  the  subspace  for  SIMO  models  is  well  established,  we  focus  on  the  novel  method: 
the  root  selection  method. 

To  test  the  performance  of  the  proposed  algorithm,  a  suite  of  simulations  were  run, 
with  the  following  set  of  parameters,  unless  otherwise  stated:  total  #  sub-carriers,  N  =  64; 
#  virtual  sub-carriers  (NSC’s)  15;  SNR  =  15  dB;  3-ray  channel  with  independent  Rayleigh 
fades,  exponential  power  delay  profile  with  decay  parameter  P  =  1/5.;  2  receive  antennas; 
M  =  10  OFDM  blocks  for  channel  estimation.  Bias  and  MSE  were  estimated  from  a 
suite  of  500  Monte  Carlo  runs;  channel  parameters  were  drawn  independently  for  each 
run.  Mean  bias  and  MSE,  for  the  set  of  L  -1-  1  channel  taps  are  reported;  in  order  to 
fix  the  scalar  gain  ambiguity  -  present  with  any  blind  channel  estimation  method  -  we 
normalized  the  channel  response  to  unit  norm.  The  symbol  sets  tested  were  16-QAM, 
64-QAM,  QPSK,  8-PSK,  16-PSK,  and  64-PSK.  Solid  lines  correspond  to  the  PSK  symbol 
sets,  and  dashed  lines  to  the  QAM  symbol  sets. 

Figure  2.1  shows  performance  vs.  SNR;  notice  that  QAM  performance  has  the  usual 
floor  effect;  PSK  performance  does  not  suffer  from  this,  and  PSK  performance  does  not 
depend  upon  the  size  of  the  alphabet.  Figure  2.2  shows  performance  as  the  number  of 
OFDM  symbols  was  increased;  with  PSK  symbol  sets  good  performance  is  achieved  even 
with  5-10  OFDM  symbols.  Figure  2.3  shows  results  as  the  actual  channel  length  was 
increased  (with  N  =  64,  and  15  NSC’s;  adequate  CP  was  used  in  every  case).  Note 
that  degradation  in  performance  is  graceful.  Next,  keeping  L  =  3  and  N  =  M  fixed, 
the  number  of  band-edge  NSC’s  was  increased;  figure  2.4  indicates  that  degradation  is 
graceful. 


2.8  Discussion 

From  the  simulations,  we  verify  that  performance  for  PSK  symbol  sets  does  not  suffer 
from  the  floor  effect,  and  that  good  channel  estimates  are  obtained  with  as  few  as  5  to  10 
OFDM  symbols. 

In  the  absence  of  CP,  spatial  diversity  has  been  exploited  for  estimating  channel  pa¬ 
rameters;  comparing  with  results  in  [38],  we  note  that  the  absence  of  CP  cannot  be 
easily  made  up  with  spatial  diversity.  Spatial  diversity  (here,  the  inclusion  of  one  extra 
receive  antenna)  tremendously  decreases  the  number  of  OFDM  symbols  required  for  re¬ 
liable  channel  estimation  by  SISO  approaches  (e.g.,  the  cyclic  approach  of  [28]  and  the 
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finite-alphabet  approach  of  [32]  both  require  substantially  more  OFDM  symbols). 

The  performance  of  the  estimators  proposed  in  Section  6  can  be  readily  accomplished 
using  the  results  of  Section  5,  along  the  lines  of  [36]. 
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Ntot=64  Nvsc=  *  snr=20  dB  L=3  M=10 16,64QAM  4,8,16, 64PSK 


#NSCs 

Figure  2.4:  Bias  and  MSE  vs.  #  null  sub-carriers 
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